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Abstract. Global entropy solutions in BV for a scalar nonlocal conservation 
law with fading memory are constructed as limits of vanishing viscosity ap- 
proximate solutions. The uniqueness and stability of entropy solutions in BV 
are established, which also yield the existence of entropy solutions in L°° while 
the initial data is only in L°°. Moreover, if the memory kernel depends on a 
relaxation parameter e > and tends to a delta measure weakly as measures 
when e — » 0+, then the global entropy solution sequence in BV converges to 
an admissible solution in BV for the corresponding local conservation law. 



1. Introduction and Main Theorems 

We study global entropy solutions to a scalar nonlocal conservation law with 
fading memory: 

u t + f(u) x + f k(t-T)f(u(T)) x dT = Q, 168, (1.1) 

Jo 

and initial data 

u(0,a)=«o(a;), (1.2) 

where / : R — > K is a smooth function and uq 6 BV(M). For simplicity, we 
sometimes use the notation u(t) :— u{t,x) to emphasize the state at t > as in 

In one-dimensional viscoelasticity, hyperbolic conservation laws 

U t + F(U) x = (1.3) 

correspond to the constitutive relations of an elastic medium when the value of the 
flux function F at (t,x) is solely determined by the value of U(t,x). However, this 
model HO|l is inadequate when viscosity and relaxation phenomena are present. In 
that case, the flux function depends also on the past history of the material, i.e. on 
U(t, x) for t < t. Under these circumstances, we say that the material has memory. 
An important class of media of this type are materials with fading memory, which 
correspond to the constitutive relations with flux functions of the form 

F(U(t,x))+ f k(t-r)G(U(r,x))dT, (1.4) 
Jo 

where F, G are smooth functions and A; is a smooth kernel, integrable over R + := 
[0, oo). When the kernel satisfies appropriate conditions motivated by physical 
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considerations, the influence of the memory term in (|1.4J) reflects a damping ef- 
fect. Consequently, global smooth solutions exist for given small initial data (cf. 
Renardy-Hrusa-Nohel |15|). in contrast to the situation with elastic media in which 
classical solutions in general break down in finite time even when the initial data 
is small. However, when the initial data is large, the destabilizing action of nonlin- 
earity of the flux function / prevails over the damping, and solutions break down 
in a finite time; see Dafermos j^J and Malek-Madani-Nohel 10 . 

In this paper we first construct global entropy solutions in BV to the nonlocal 
conservation law (|l.l|l with fading memory via the vanishing viscosity approxima- 
tion by adding the artificial viscosity term as follows: 

< + /(«")*+/ k{t-r)f{u' J {T)) x dT = v{u» + f k{t-T)u v {r)dr) xx . (1.5) 
Jo Jo 

For a scalar local conservation law modeling elastic materials, such a result was 
first established in P OUT] . 

The main motivation for the vanishing viscosity approximation (|1.5|) is that the 
conservation law (|l.l(l can be viewed as a linear Volterra equation, which was first 
observed by MacCamy ^1] and later employed in Dafermos |1] and Nohel-Rogers- 
Tzavaras (also see pQ). In this way, it is easy to extract the damping character 
of the memory term. Let r(t) be the resolvent kernel associated with k(t): 

r + k * r = —k. (1-6) 

Then we can write (|1 . 1|> as 

-f{u)x=u t + r(t - t) v! (t) dr. 
Jo 

Integrating by parts yields 

u t + f(u) x +r(0)u = r(t)u Q - f r' (t - r)u(r) dr, (1.7) 

Jo 

which is equivalent to (|l.l(l . The vanishing viscosity approximation 1)1.5(1 is equiv- 
alent to the following artificial viscosity approximation to (|1.7|) : 

u v t +f[u v ) x + r{<S)u v =r{t)u Q -f r' (t - t)u v (t) dr + v u v xx . (1.8) 

Jo 

Note that the above argument applies only if the nonlinearity / in the instantaneous 
response is the same as in the memory term chosen in (|1.1|) ; also see 0] El Q2| f° r 
the same restriction. The artificial viscosity term in (|1.5fl is chosen so that (|1.8|) 
has the standard artificial viscosity term to ensure the L°° and BV estimates of u v 
(also see 00] H2]). 

The existence of a unique, regular local solution u v (t, x) of (|1.8|) when the initial 
data uo is smooth can be established through the standard Banach Fixed Point 
Theorem. The local solution may be extended to a global solution with the help of 
the apriori L°° estimate established in Section 2.1. 

A function u — u(t, x) is called an entropy solution to the Cauchy problem 
(HHl-JOl if it satisfies that, for any test function ip £ C„ (K+) with R 2 + := R+ x R, 
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T - > 0. 

(j](u)(pt + q(u)ip x + r]' (u)(r(t)uo — r(0)u — / r (t — T)u(T)dT)ip) dtdx 

Jo 

r){u (x))ip(0,x)dx > 0, (1.9) 

for any convex entropy rj(u), where q(u) is the corresponding entropy flux satisfying 
q'(u) = i{u)f{u). 

Before we state the results, we introduce some notations. Let g be the standard 
mollifier. We define the mollification of uo to be 

Uq := (uqXv) * Qu, (1-10) 

where, for each v > 0, Q v ix) :— i g(^) and Xv( x ) := 1 f° r M < 1/^ an d otherwise. 
The main result is the following. 

Theorem 1.1. Consider the Cauchy problem (|1.5fl with Cauchy data: 

u' / (0,x)=uZ(x), (1.11) 

where the initial data Uq is given by <|1.1(J|> anduo 6 BV(M). Let the resolvent kernel 
r associated with k as defined in l|1.6f) be a nonnegative, non-increasing function in 
L 1 (R + ). Then, for each v > 0, the Cauchy problem (|1.5f) and (11.111) has a unique 
solution u v defined globally with a uniform BV bound. Moreover, as v — > 0, u v 
converges in Lj oc to an entropy solution u £ BV to (ll.l|) - <|1.2f) , which satisfies 

NIl°°(r*) < II w o||l~(r), (1-12) 
TV{u(t)}+ [ r(t-T)TV{u(r)}dT<CL, (1.13) 



||«(t)-«(*)||£l(R) < C\t-s\, (1.14) 

where L = 1 + M|z,i(r + ), and C = C(TV{uq}, \\uo\\l°°) is a positive constant 
independent ofr(t). 

Furthermore, we have 

Theorem 1.2. Let the resolvent kernel r(t) associated with k be a nonnegative and 
nonincreasing function mi 1 (IR+). Letu,v G BV(M.\) be entropy solutions to (|1.1|) 
with initial data uq,vq £ BV(M.), respectively. Then 

\\u(t) - v(t)\\ L i m + / r{t - t)\\u(t) - v(t)\\ l i (m) dr < L \\u - v \\ L i (M) . (1.15) 
Jo 

That is, any entropy solution in BV to (ll.l|) - (|1.2H is unique and stable in L . As 
a consequence, if Uq is only in L°° , not necessarily in BV(M), there exists a global 
entropy solution u £ L°° to (|l.lfl - l|1.2l) . 

Having established the above results, we then analyze the case when the kernel in 
the scalar equation (|1.1|) is a relaxation kernel k £ that depends on a small parameter 
e > so that k £ (t) — i (a — 1) 6(t) weakly as measures when e — > 0+, where 5(t) 
denotes the Dirac mass centered at the origin. That is, 

U £ t +f(u e ) x + f k £ {t-T)f(u e (T)) x dT = 0, (1.16) 



m e (0, x) = u Q (x) 6 BV(R) (1.17) 
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with sup e>0 ||fe e ||ii( R+ ) < oo. 

We denote the entropy solution to the above problem by u £ (t,x). Let r £ be the 
resolvent kernel associated with k £ via (|1.6I) . Hence, l|1.16|) reduces to 

u E t +f(u £ ) x + r £ (0)u E = r £ (t)u - f r' £ (t - t)u(t) dr (1.18) 

Jo 

with sup e>0 Hrellii^) < oo. 

By Theorems ll.llfTT^l we conclude that the unique entropy solution sequence 
u £ G BV to l|1.16|l ~ H1.17|l is uniformly bounded in L°° and uniformly L 1 -stable, 
independent of e. Then the solution sequence {u e } is a compact set in Lj oc so that 
we can extract a subsequence {u 6k } that converges in L\ oc to an admissible weak 
solution of the local conservation law: 

ut + af(u) x = 0, (1.19) 

with Cauchy data u £ BV. 

Theorem 1.3. Consider the Cauchy problem <|1.16|) - <|1.17fl with uq £ BV(M). Let 
the resolvent kernel r £ associated with k £ as defined in (|l.tj|) be a nonnegative, 
nonincreasing function with uniform L x -norm independent of e. Then the entropy 
solutions u e to (|1.16fl - <|1.17|) are uniformly bounded in L°° and stable in L : 

\\u e (t)\\ L ^ {M 2 +] < ||u || L oo (K) , 

TV{u £ {t)} <CL, 

\\u £ (t)-u'(s)\\ L i {R) <C\t~ S \, 

\\u e (t) - v e (t)\\ L i {m < L\\u - uo||li(r), 

where L := 1 + sup e>0 < oo, C = C(TV{uq}, \\uo\\l<*>) > is independent 

of e, and v e (t, x) is the entropy solution to l|1.16(l (| 1 . 1 T|l with initial data vq 6 BV . 
Furthermore, if k £ (t) — 1 (a — 1) <5 (t) weakly as measures when e — > 0+ ; then u e 
converges in Lj oc to an admissible weak solution u of the Cauchy problem (|1.19|) 
and p.2|l with initial data uq G BV. 

In Theorems II . II fOl the assumptions on r £ (t), or r(t), can easily be converted 
to the assumptions on k £ (t), or k(t), because of their symmetry between the kernel 
and the resolvent through ill.6|) . For example, such kernels k £ (t) especially include 
the following Set of Kernels (i)-(vi): 

(i) k' £ (t) > and ||fc e ||Li(E + ) < K for some constant K independent of e > 0; 

(ii) |1 + k £ (z)\ ^ for any z with Re(z) > 0, and k £ (z)(\ + k £ (z)) < for the 
Laplace transform k £ of k £ ; 

(iii) sup ue R|(l + k £ (iuj))~ 1 \ < q for some constant q independent of e; 
(vi) There exist positive numbers T ~ e and t ~ e such that 

/ M*)l<4-. su p I \k s {t)-k £ {t-s)\dt<-. 

J\s\>T LZ Q 0<s<tiJr 4 

The prototype is 

1 — a t 

k B (t) = exp{ — ), < a < 1, (1.20) 

e e 

for which the corresponding family of resolvent kernels is 

r £ (t) = l -^ eX p(--). (1.21) 
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In Section 2, we develop techniques for the nonlocal case, motivated by VoPpert- 
Kruzkov's techniques |17l l%] and the L°°-estimate techniques for the local case, to 
establish uniform L°° and BV estimates of the vanishing viscosity approximate 
solutions, independent of v, by using the damping nature of the memory term. 

As a corollary of these estimates, we establish the convergence of the vanishing 
viscosity approximate solutions to obtain the existence of entropy solutions in BV. 
In Section 3, we show that the entropy solution in BV is unique and stable in L 1 
with respect to the initial perturbation. In Section 4, we prove Theorem 11.31 and 
discuss the hypotheses of the theorems. Finally we give an example and show the 
relation of the fading memory limit with the zero relaxation limit as first considered 
systematically in Chen-Levermore-Liu 0; also see [HlEl^l for the model. 

2. Proof of Theorem 11.11 

In this section, we establish the uniform L°° and BV estimates, as well as the 
uniformly continuous dependence on time in L 1 , which are used not only for the 
global existence of the vanishing viscosity approximate solutions, but also for their 
compactness. We also establish the existence and regularity of entropy solutions. 

2.1. L°° Estimate. We first obtain a uniform L°° estimate. Note that, by em- 
ploying the resolvent kernel r, equation (|1.5|) can be written in the form (|1.8(l . i.e. 
we study the Cauchy problem 

< + /K) a +r(0)u"=r(tK- f 1 r\t-T)u"{T)dT + vu v xx , (2.1) 

Jo 

u(0,x) = u%(x), xeR, (2.2) 

where Uq(x) is defined in (|1.1Q(I . By rescaling the coordinates, (t, x) — > {s,y) — 
(rot,rax), we rewrite ljl.8|) as 

— r \ — ) u 2 / r ( 

where ro > r(0) is any positive constant and u u (s,y) := u v {^-, ^-). For any even 
integer p, multiplying l|2.3|l by p \u I/ \ p ~ 1 and integrating over [0, S] x R, we obtain 

/>S />oo 

\u"(S)\"dy+p / \&(8,v)\*dyds 

JO J-oo 

/oo r-S roo -, 

\<(v)\ p dy+p / -r{-)\ul{y)\\u»{s,y)\P- x dyds 
~oo JO J-oo r o r 

-P / |«"(a,J/)r 1 / -^r'(—W(r,y)\dTdyds. (2.4) 

JO J-oo Jo r r 

By employing the standard inequality ab < eo + £i y with £o = ( 4( ' :p p ~ 1 ' > 

e i = 4(p-i) ano - p + q = 1' ^e second term on the right-hand side of Ij2.4|l is 

estimated as 

P / -K— ) \<iv)\ \* V (*,V)\ P dyds 
Jo J-oo r o r a 

nS f OO / 1 „ \ P 1 r ^ r °° 



u: + .f(n u )y+^ = ~r (-)<-- r'( )«"(t) dr + ur u v yy , (2.3) 



10 J-oo 



< / eo -r(-) \u»(y)\P dyds + -p / \u v (s,y)\* dyds, (2.5) 



ro r J 4 
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so that the last term in Ij2.5|l is dominated by the damping in Ij2.4|l . Similarly, we 
treat the last term in (|2.4J) to obtain 

\u v (S,y)\*dy< r H(y)\ p dy + e f* ( -r(-)Y ds f°° \u» (y)\ p dy 

J-oo Jo \ r Q r J J-oo 

+ P(S) / / \u v ( T ,y)\ p dydT, 



q \p/q 

dr ds. By Gronwall's inequality, we have 



where /3(5) := e J Q S (/„' | ^r'(^) 

rS poo pS 

/ \u v (r)\ p dydT <e$° ^ dw / W(s)ds (2.6) 

for W( S ) := (l + eo/ * ds) J \u» (y)\v dy. Let 

if s = SU p(l|r'(-)|: 5€[0,5]}. (2.7) 
Then, for all s E [0, S), we have 

P(S)<e K s [-- / r'( )dr < e K s 5, (2.8) 

1 /- s 

since r is nonincreasing, and hence lim — / (3{w) dw = 0. Also, 

p^oo p J 

(f wt - )d -) " s si (' +4 (f (H>)'*) ') ' ■ 

Thus, if we raise l|2.t)[) to 1/p, take the limit as p — > oo, and note that £o(p)' ? - > 4 
as p -> oo, we conclude that ||u I/ || L oo([ ^ ]xK ) < ||(1 + ^r(^))|tZo(j/)||U«([o,S]xR)- 
That is, 

||«"||£-(r») < (1 + -||r|| i o 0(K p)|HU»( K) . (2.9) 

Since r(t) > is nonincreasing so that < r(t) < r(0) < oo and ro > r(0) is an 
arbitrary constant, then, as ro — > oo, we conclude |j?2 ly || i oc( R 2 j < ||wo||l oc (r)- By 
rescaling the coordinates backwards, we obtain the uniform L°°-bound on u" to 
(|1.8|l independent of v; 

IKIIl°°(r*) < ll"olk-(R)- (2-10) 

2.2. BV— Regularity and Estimates. First note that Uq £ C°° has the following 
bounds: 

||(«o)x|Ux < M(« ), IKOxx]]^ < ^Af(tio), (2.11) 

for some Co > independent of f, where M(u ) := TV{uq} + 2||mo||l° c - 
Set v = u v x and u> = u\. Then the evolution equations of v and w are 

1H + (/'(«*>)« + r(0)v = r(t)«) x - [ r'(t - t)v(t) dr + vv xx , (2.12) 



v(0,x) = (u v ) x , (2.13) 
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and 



w t + (f'{u u )w) x + r(0)w = - { r'(t - t)w(t) dr + vw xx , (2.14) 



w(0, x) = <(0, a;) = «) M - /'«) (u v ) x . (2.15) 
Multiplying (|2.12|) by sgn(v(t,x)) and integrating with respect to x, we get 



d 



M n v(t)\\ L i + r(0)||«(t)|Ui < r(t)M(u ) - J r'(t - t)\\v(t)\\v dr 
since r(0) > and r is nonincreasing. Integrating over t 6 [0, T] yields 

\\v(T)\\ L1 +r(0) I* \\v(t)\\ L idt 



<(l+/ r(t)dt)M(u ) - [ ( r'(t-T)\\v{T)\\ L idTdt. (2.16) 
Jo Jo Jo 

Changing the order of integration in the last term, we arrive at 

\\v(T)\\ Ll +r(0) f \\v{t)\\ L ,dt 
Jo 

<(1+/ r(t)dt)M(u )- [ (r(T-r)-r(0))\\v(r)\\ L1 dT. 
Jo Jo 

Thus, we have 

\\v(T)\\ L i+ [ r{T - t)\\v{t)\\ l , dr < (1 + f r(r) dr)M(u Q ). 
Jo Jo 

Because r(-) is bounded in L 1 (M + ), we obtain the following uniform bound on the 
gradient v = u x , 

\\<{t)\W+ f 'r(t-r)\\u^r)\\ L rdT<LM(u ), (2.17) 
Jo 

where L := 1 + |M|li(r + ). Similarly, we have 

IK*)IUi+ [ r(t-T)\\w(T)\\ L .dT<\\w(0)\\ L i. (2.18) 
Jo 

Using (|2.1|1 and the bounds in (|2.11() for the initial data, we find from 12.15fl that 

IM0)|| L i < v\\(u v Q ) xx \\ L i + \\f(u v Q )\\ Laa M(u ) < CiM(«o) 

for some C\ > independent of v. Hence, by l|2.18|l . is uniformly bounded in 
L . Thus, for < s < t, we get the uniformly continuous dependence on time for 
the solutions to l|2.1 |l -l|2.2 |l : 

Htt^t) - u"(a)||ii < J \\w(r)\\ dr <C\t-s\ (2.19) 
with C = max(Ci, 1) M(uq) independent of v. 



8 



GUI-QIANG CHEN AND CLEOPATRA CHRISTOFOROU 



2.3. Existence of Entropy Solutions in BV to ifO jl -lfO jl ■ Using (JOB, f2~T7jl . 
and (|2.19(l . Helly's Compactness Theorem yields that a convergent subsequence 
{u"'"} may be extracted with v m \ as m — > oo, whose limit is denoted by u, i.e., 

u Vm {t) — >u(t) mLj oc forallt>0. (2.20) 

The limit u{t, ■) is a BV function satisfying ((i"T2l - l(i"Xl|) for all t, s > 0. By 
construction, it is easy to check that the limit function u(t, x) is an entropy solution 
to fl~Tl)-(01). 

3. Proof of Theorem 1.2 

In this section, we prove the uniqueness of entropy solutions in BV as stated in 
Theorem ll.2l For any u £ BV, the whole space can be decomposed into three 
parts (see 0E1[T7]): 

R 2 + = J(u)UC{u)Ul{u), 

where J(u) is the set of points of approximate jump discontinuity, C{u) the set of 
points of approximate continuity of u, and I(u) is the set of irregular points of u 
whose one-dimensional Hausdorff measure is zero. 

First, the entropy inequality (|1.9|l implies that, on a shock x = x(t) in J(u), 

a[ V (u)} - [q(u)} > 0, (3.1) 

where [ri(u)] = rj(u(t, x(t) + 0)) — r](u(t, x(t) — 0)) and a = x'{t) is the shock speed. 

Now assume that u,v 6 BV(M.+ ) are the entropy solutions with initial data 
Uq, Vq £ BV(R), respectively. Then it can be easily checked that, on J{u) U J(v), 

a[\u - v\] - [sign(u - v)(f(u) - f(v))] < 0. 

In the continuous region C(u) P\C(v), since r'(t) < 0, 

H(t,x) := \u(t)-v(t)\ t +q(u,v) x +r(0)\u(t)-v(t)\ 

- r(t)\u - v \ + f r'(t-T)\u(r) -v{r)\dT 
Jo 

t 

\r'(t - t)| (\u(t) - v(t)\ - sgn(u(t) - v(t))(u(r) - v(t))) dr < 0. 
Therefore, [i as a measure on satisfies 

»(R 2 + ) = - {<rM-[q]) + tiC(u)nc(v))<o. 

J{u)UJ(v) 

Then we follow the same steps as for the BV estimates in Section 2.2 to conclude 
(I1.15fl . When uq S L°°, let Uq be a sequence of initial data in BV for which Uq — » uq 
as k —> oo. Then the ^-stability result (|1.15(l implies that the corresponding 
entropy solution sequence u k £ BV to Ijl.lfl with data Uq is a Cauchy sequence in 
L 1 which yields a subsequence converging to u(t, x) £ L°°. It is easy to check that 
the limit u(t, x) is an entropy solution. 
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4. Proof of Theorem 11.31 

Let u s £ BV denote the unique entropy solution to 1)1. 18[) with initial data 
uq £ BV. Then the solution sequence {u 6 } is uniformly bounded and is uniformly 
stable in L 1 with respect to the initial data since 

IKIL-(R 2 + ) < IKIU°°(R), \\u £ (t) - V £ (t)\\ L i < L \\u - V \\ L i 

and satisfies the following apriori uniform bounds: 

TV{u e (t)}<CL, \\u e (t)-u e (s)\\ LHm <C\t-s\. (4.1) 

This implies that there exists a convergent subsequence {u Em } with e m — * as 
m — > oo, whose limit is denoted by u, i.e., u £m (t, x) — ► u(t,x) in L; oc . Then, since 
fc s (t) — 1 (a — l)5(t) weakly as measures when e — > 0, we conclude that u is an 
admissible weak solution of the Cauchy problem (|1.19|) and (|1.2fl . The proof of 
Theorem II .31 is complete. 

Finally we discuss some families of kernels {k E } that satisfy the assumptions 
stated in Theorem 1 1.31 

Suppose that k £ £ L 1 (Ht+) for all e > 0. Then, by the Paley- Wiener Theorem 
|14j. the resolvent r e of k £ is in L 1 (M + ) if and only if 

\l + k s (z)\^Q for all Re(z) > 

for the Laplace transform k s of /c e . By extending k £ as zero to the negative real 
axis, we choose the number q such that 

q > snp weK |(l + k e (iu))~ 1 \, 

and choose positive numbers T and 77 satisfying 

/ %{t)\dt< Sup / |* e (*)-*e(*-*)l*<7- 

J\s\>T 1Z Q 0<s<r]J-oo 4 

Then 

\\r s \\ L HM + ) < (8r6gT||fc £ || il(R+) l \8\\k £ \\v (M+) /r,] + 6) q\\k £ \\ LHR+h (4.2) 

where [s] denotes the smallest integer > s. 

With this, for each e, we can take k e £ L 1 such that |1 + k e (z)\ ^ for all 
Re(z) > and the numbers defined above to be: q independent of e, T ~ e and 
77 ~ e. Then, by l|4.2|l . r e is uniformly bounded in L 1 (R). Furthermore, any kernel 
fc £ (t) in the Set of Kernels (i)-(vi) satisfies the assumptions in Theorem 1.2. 

A prototype is the family of kernels k £ (t) in (|1.20fl that satisfies these assump- 
tions. Then the corresponding family of resolvent kernels is r e (t) in p. 21(1 which 
fulfills the assumptions of Theorem II .31 when < a < 1. For this example, the 
scalar nonlocal equation (|l.l|l : 

1 — a r f t 

u t + f(u) x / e ~J(u(t)) x dr = 

e Jo 

can also be written as a system of two equations: 

ut + (f(u) - v) x = 0, 

(l-a)f(u)-v (4.3) 
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Then the range of a £ (0, 1) is the sub-characteristic condition. Thus, the re- 
sult of Theorem 11.31 applying to this special case is equivalent to establishing the 
convergence of the relaxation limit l|4.3[l as considered in [31 ^3 ^] . 

Remark 4.1. In order to obtain that the resolvent r e of k e is integrable and 
ll r e||L 1 (R+) ^ 20 7 it also suffices by the Shea-Wainger Theorem \7\ to require that a 
family of kernels {k £ } satisfies that, for each e > 0, k £ £ Lj oc (R + ) and is nonneg- 
ative, nonincreasing, and convex on R + . 
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